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Zusammenfassung — £s wird cin Verfahrcn zur hochaufldsen- 
den, weUenformerhaltcnde Digltalisierung analogcr Slgnale vor- 
gestellt Die fibliche skalarc logarithmische Quantisierung wird 
mehrdimensionale Kugelkoordlnaten tibertragen und die 
us resultierenden Vortcile wie z.B. cin konstanter Stftrab- 
Qber cinen cxtrem hohen Dynamikbcreich werdcn dls- 
iert Um die im Quellensignal vorhandcnen statistlschen 
hhanglgkeitcn fQr eine weitcren Stdrabstandsgcwinn ausnutzen 
zu kttnnen, wird einc Methode vorgcstellt, Differentielle Puisco- 
demodulation (DPCM) mlt spharisch logarithmischer Quantisie- 
rung zu kombinicren. Das resuitierende Verfahren erzlelt eine 
wirksame Datenreduktion bei extrem geringcr Signalverzoge- 
rung. 

I. ElNLBITUNO 

Viele Bereiche der digitalen Nachrichtenubertragung haben 
das gemeinsame Ziel, ein zunachst analog vorliegendes Quel- 
lensignal (MeBwerterfassung, Audiosignalverarbeitung in Stu- 
dioqualitat, usw.) unter Ausnutzung der Vorteile der digitalen 
Obertragung zu einem EmpfSnger zu transportieren, um es 
dort in Form eines analogen Ausgangssignals wiederzugeben. 
HierfUr muss das Analogsignal zum einen digitalisiert bzw. 
quantisiert werden, zum anderen muss eine entsprechende 
Datenkompression erfolgen, um das Signal bei mOglichst guter 
gualitat Uber ein digitales Obertragungssystem mit begrenzter 
zu tibertragen. Die zahlreichen bereits existierenden Ver- 
bn zur Digitalisierung analoger Quellensignale verfolgen 
Wise grundsatzlich verschiedene Ansatze, sie lassen sich 
jedoch alle in zwei Hauptkategorien einteilen: wellenformer- 
haltende und nicht wellenformerhaltende Verfahren. Letzte- 
re sind beispielsweise im Audiobereich von entscheidender 
Bedeutung, da hier oftmals nur der subjektive Horeindruck 
am Empfangerausgang entscheidend ist (z.B. Ausnutzung psy- 
choakustischer Maskierungseffekte). Hier spielen Wellenform- 
verfinderungen (Amplituden- und Phasenverzerrungen), die 
das Signal durch Quantisierung und Komprimierung erfShrt 
i.d.R. nur eine untergeordnete Rolle. Die Obertragungsqualitat 
eines solchen Verfahrens lSsst sich nicht durch einen St5rab- 
stand 101og 10 (SNR), SNR: Signal to Noise Ratio, im klas- 
sischen Sinne messen; sie muss vielmehr durch aufwendige 
HOrtests mit trainierten Personen ermittelt werden. Fttr viele 
Anwendungsbereiche (z.B. in der Messtechnik) sind derartige 
Signalcodierungsverfahren jedoch ganzlich unbrauchbar. Diese 
Arbeit befasst sich mit den wellenformerhaltenden Verfah- 
ren, deren Hauptziel eine mBgiichst unveranderte Wellen- 
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form am Decoderausgang und somit ein m&glichst hotter 
objektiv messbarer Stttrabstand ist, vorzugsweise im Sinne 
des Verhaitnisses von Varianz des Nutzsignals zum mittleren 
quadratischen Fehler. Hier wird ein Verfahren vorgestellt, das 
zudem Uber eine extrem hohe Dynamik verfUgt (d.h. das SNR 
ist konstant Uber einen sehr groBen Aussteuerbereich von z.B. 
60 dB und mehr) und bei nur extrem geringer VerzBgerung 
in der Lage ist, sowohl statistische Abhfingigkeiten in der 
Folge der Abtastwerte aus dem Quellensignal als auch die 
gUnstigen Packungseigenschaften mehrdimensionaler Gitter 
(Vektorquantisierung) zu nutzen. 

II. LOOARITHMISCHB QUANTISIERUNG 

Es wird eine hohe Dynamik eines Quantisierers, also ein 
weiter Bereich des mittleren Signalpegels angestrebt, in wel- 
chem das SNR und somit der maximale relative Quantisie- 
rungsfehler konstant ist Hierbei bezeichnen die 

Breiten der Quantisierungsintervalle und r< die jeweiligen 
Rekonstruktionswerte. Diese Aufgabe fUhrt bekanntermaBen 
zu der logarithmischen Quantisierung zJB. im Sinne des 
.A-Laws, wie sie in der Fernsprechtechnik schon seit langem 
angewendet wird. FUr den Bereich mittlerer Aussteuerung 
ergibt sich fur R » 1 folgender St6rabstand [3]: 



101og 10 (SNR) = R- 6,02dB + 101og 10 
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(1) 



wobei R die mittlere Rate der Quantisierung 
(bit/Abtastwert) bezeichnet und A den Ublichen Parameter der 
logarithmischen Quantisierung darstellt; d.h. der maximale 
relative Fehler ist fUr Abtastwerte mit einem Betrag > l/A 
(bezogen auf einen Signalbereich von -1 bis +1) konstant. 
Damit wird ab einem mittleren Signalpegel des Quellensignals 



von 



-Bi := 201og 10 (lM)dB 



(2) 



ein StSrabstand gemaB (1) erreicht. Der Wert Bi charakte- 
risiert damit den Dynamikbereich der logarithmischen Quan- 
tisierung. Innerhalb dieses Dynamikbereiches ist der Storab- 
stand zudem ganzlich unabhangig von der Wahrscheinlich- 
keitsdichtefunktion (WDF) des Quellensignals so dass das 
Verfahren universell einsetzbar ist Allerdings steilt der Term 



101og 10 (3/(1 + ln(.4)) a ) im Vergleich zu ciner gleichmttfii- 
gcn Quantisierung bei Vollaussteuerung durch ein gleichvcr- 
teiltes Signal gewissermaBen einen "Verlust" durch Kompan- 
dierung dar, der filr einen groBen Dynamikbereich zu bezahlen 
ist 

III. QUANTISIERUNG IN KUOBLKOORDINATBN 
A. Kugelkoordinaten in D Dimensionen 

Aus dem Gebiet der Vcktorquantisierung ist bekannt, dass 
durch eine kugeldichte Packung in vielen Dimensionen be- 
trUchtlichc SWrabstandsgewinne bei der Quantisierung erzielt 
werden, auch wenn keine statistischen Abhfingigkeiten inner- 
halb der Signale ausgentitzt werden kttnnen. Bei sphMrischer 
Quantisierung werden D Abtastwerte (x := Xi.xa, ,.;»xi?) 
zusammengefasst, urn einen Punkt im JR D mittels Polarkoor- 
dinaten zu beschreiben. Der Radius r sowie die D - 1 Winkel 
durch folgende Gleichungen gegeben (j: imagin&re Ein- 




r 2 = 



<Pi = 



D 

1=1 

arg(xi + jx 2 ) 



<Pi = arg 



€ [-7T,+7T) 



(3) 
(4) 



(5) 



Die Rekonstruktion der 
81,24, 2d erfolgt hieraus gemftB 

xi = r • 6<-i • sin 



*6{2,...,D-1}. 
kartesischen Komponenten 



% € {2,...,D} 



xi = r • 61 • cos {(pi) = r • bo 



(6) 
(7) 




wobei fQr die Radien der •'Breitenkreise" der Kugel mit 
Radius 1 gilt: 

&z>-i = 1, (8) 

D-l 

bi = 1 • II cos (w), * € {0, D - 2} . (9) 
i=i+i 

Damit ist ausgehend von vd-i eine aufwandsgUnstige 
rekursive Darstellung gefunden. Eine aufwandsgUnstige Im- 
plementierung von Vor- und Rticktransformation ist mit Hilfe 
des sog. CORDIC-Algorithmus mbglich. 

B. Spharisch logarithmische Quantisierung 

Damit bei der Quantisierung in Polarkoordinaten die Vor- 
teile der UnabhSngigkeit des SNR von der WDF des Quel- 
lensignals weiterhin genutzt werden kflnnen, wird nun der 
Betrag (Radius) logarithmisch gemaB den Vorschriften des 
.A-Law quantisiert Trotz einer gleichm&Bigen Quantisierung 
der D - 1 Winkel (z.B. im 3-dimensionalen Raum von 
Azimuth- und Elevationswinkel) sind aufgrund der Proportio- 
nality der LSnge eines Kreisbogensegmentes zum Radius die 
Anforderungen an die logarithmische Quantisierung damit in 
alien D Dimensionen erftillt, obwohl nur in einer Dimension 



(bei der Quantisierung des Betrags) der bereits angesprochene 
Verlust durch logarithmische Kompandierung in Kauf genom- 
men werden muss. 

Urn die Implementierung der Quantisierung in Kugelko- 
ordinaten zu vereinfachen gehen wir von nttherungsweise 
kubischen Quantisierungszellen aus. d.h. die Oberflfiche einer 
P-dimensionalen Kugel mit Radius 1 wird durch ein Netz 
D - 1-dimensionaler (Hyper-)WUrfel quantisiert. Auf diese 
Weise wird eine sehr einfache Implementierbarkeit von Quan- 
tisierung und Signalrekonstruktion erreicht, wobei durch die 
iterative Vorgehensweise gem&B (4) bis (9) eine schrittweise 
Verarbeitung der einzelnen Koordinaten wie im Fall skalarer 
Quantisierung bewahrt wird. Wie nachfolgend gezeigt wird, 
ergibt sich durch diese suboptimale Quantisierung der Kuge- 
loberflache im Vergleich zu einem optimal dicht gepackten 
Gitter jedoch im Grenzwerl D 00 lediglich ein Verlust 
von 101og 10 (x) 53 1*53 dB, also ein Ratenverlust von ca. 
l/4bit/Abtastwert, der zugunsten einfacher Implementierbar- 
keit hingenommen wird. Asymptotisch wird damit die Grenze 
der Rate-Distortion-Funktion fttr die unabh&ngig und identisch 
verteilte (i.i.d.) gauB'sche Zufallsvariable bei zugleich gegen 
unendlich strebendem Dynamikbereich fUr konstanten StCrab- 
stand nur urn 1/4 bit bzw. 1,53 dB verfehlt. 

Da im Sinne der Zielsetzung einer logarithmischen Kom- 
pandierung bewusst von keiner Kenntnis einer Wahrscheinlich- 
keitsdichtefunktion fUr i.i.d Quellensignalwerte ausgegangen 
wird, stellt aber wiederum die Rate-Distortion-Funktion fUr 
die gauB'sche Zufallsvariable die untere Schranke fUr die 
erreichbare Rate (bzw. Verzerrung) bei gegebener Verzerrung 
(Rate) dar (Berger's upper bound der Rate-Distortion-Function 
[1]). Somit sind unter den hier gestellten Anforderungen 
und EinschrSnkungen ttber 1,53 dB bzw. 1/4 bit/Abtastwert 
hinausgehende Verbesserungen ohnehin nicht mehr mSglich. 

Die Aufteilung der M D (d.h. M := 2 R ) pro Quantisie- 
rungsschritt zur Verfiigung stehenden Quantisierungsstufen auf 
die einzelnen Dimensionen (also auf den Radius und die D — 1 
Winkel) wird auf folgende Weise vorgenommen: 

FUr den logarithmischen Bereich ^ < r < r 0 der .4-Law- 
Vorschrift gelte die folgende Kompressorkennlinie 

k(r) = r 0 • (c • In y- + 1^ mit c := — : 1 

bzw. fUr deren Ableitung 



+ ln.A 



(10) 
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wobei A den freien Parameter des .A-Law bezeichnet und tq 
als noch nSher zu beschreibender Normierungsfaktor fur den 
Radius eingefUhrt wird. Gehen wir nun von M v Quantisie- 
rungsintervallen fur den Radius aus, so gilt fUr die Breite der 
Quantisierungszelle in Richtung des Radius 

Ar(r) m ro 1 



M t • Jk'(r) M t • c 



(12) 



Man beachte, dass Ar(r) im betrachteten Bereich unabhangig 
von ro ist 



Urn bci der Quanlisicrung in Polarkodrdinaten ein ein ahnli- 
ches Obersteuerungsverhalten zu crhalten wie bei kartesischen 
Koordinaten, normicren wlr den Maximalwert ro fiir den Ra- 
dius derart, dass die somit entstehende JD-dimensionale Kugel 
das gleiche Volumen besitzt wie ein D-dimensionaler WUrfel 
der KantenlMnge 2 (entsprechend einem Quantisierungsbereich 
X{ € [-1; 1] in jeder Dimension). 

Vkugcl = V\vucrfel 



1 2 D 



mit dem Volumen <xd der D-dim. Einheitskugel gem&B [4] 



(D/2)\ 



Damit gilt: 



ro = 




r 0 > 1 VI> € IN mit D > 2 



(13) 



(14) 



r m die Bestimmung des resultierenden Stdrabstandes zu 
vereinfachen, wird im Folgenden r = 1 gesetzt, d.h. es 
werden diejeriigen Quantisierungszellen betrachtet, welche auf 
der Oberflache einer Kugel mit Einheitsradius liegen. Dies 
ist ohne Beschr&ikung der Allgemeinheit mttglich, da das 
SNR aufgrund der logarithmisctien Quantisierung im Bereich 
5f < r < ro unabhSngig vom Radius ist. Fiir die Breite der 
Quantisierungszellen erh&lt man nun 

Ar(l) * — := A (15) 

M r • C 

Durch die bereits angesprochene gleichm&Bige Winkelquan- 
tisierung wird die Oberflfiche der D-dimensionalen Einkeits- 
kugel in Mk Zellen eingeteilt, welche jeweils die Form von 
D - 1-dimensionalen WUrfeln besitzen. FUr die Oberflfiche der 
Kugel gilt [41: 



woraus sich wiederum mit (IS) die Kantenlfinge der Quanti- 
sierungszellen folgendermafien darstellen lassen: 

*-*•(¥)*• ' « 

Die verbleibenden Mk Quantisierungsintervalle werden in 
Afu.i, ...Mk,D-i aufgeteilt und fOr die Quantisierung der D-l 
Winkel verwendet. Auch hier erfolgt die Aufteilung optimal 
in Bezug auf konstantes SNR: 

«... = l=#j -[«-••.• Gc)*J. 

<€{2,...,£>-l} 

mit \x\ € N: gr&Bte ganze Zahl < z mit x 6 R. 

D. Quantisierungsgerdusch und Sttfrabstand 

Fiir die Leistung des Quantisierungsgerfiusches bei D- 
dimensionalen kubischen Quantisierungszellen gilt bekannt- 
lich nfiherungsweise (M » 1) 



S = p D - r D ~ l mit p D = D • 



(16) 



in unserem speziellen Fall r = 1 folglich So = Pd- 
it A D_1 als mittlerer Fl&chenbeitrag einer Quantisie- 
gszelle zur Kugeloberflache gelte nun 



Mk • A D *" 1 = Pd 



(17) 



Somit ist ein fairer Ausgleich zwischen der Quantisierung des 
Radius und der Quantisierung der Oberflache erreicht. 

C. Aufteilung der Quantisierungsstuf en auf Radius und Winkel 

Gehen wir von M Intervallen fUr jeden der D zu quanti- 
sierenden Werte xi aus, so gilt fiir die Quantisierung in D~ 
dimensionalen Polarkoordinaten 



M k • M r = M D . 



(18) 



Mit (15), (17) und (18) erhfilt man fiir die Anzahl der fiir die 
Quantisierung des Radius zur VerfUgung stehenden Quantisie- 
rungsintervalle 



M t = M« 



C {D-1)/D 



A 2 _ 
— • D 
12 



(22) 



Aufgrund der eingangs eriauterten UnabhSngigkeit des 
Sttirabstandes vom Nutzsignalwert gilt fiir die KantenlMnge der 
Quantisierungszelle in der Kugel folgender direkt proportiona- 
ler Zusammenhang 



A(r) = A • r = Ay> ■ r 



(23) 



Somit kann fiir. die Berechnung des StGrabstandes 
(101og 10 (SNR)) ein D-dimensionaler Nutzsignalwert 
mit Leistung r 2 = 1 verwendet werden. Unter Zuhilfename 
von (20) erhfilt man 

SNR = — = F(D) . M 2 



AI 

12 



(24) 



mit 



^ (25) 



((D/2)\) 2/D -<?' D 



(26) 



und daraus mit (13) und (16) 

F(D) kann als Verlust bezQglich der Rate-Distortion- 
Schranke fUr i.i.d. gauB'sche Zufallsvariablen (6 dB-pro-Bit- 
Regel) interpretiert werden und ist in Abb. 1 fiir verschiedene 
Werte von A dargestellt. 

MitHilfederStirlingschenNSIierungx! w yfimx (f ) x ISsst 
sich der Grenzwert 



lim F(D) = — = -l,53dB 



(27) 



09) 



berechnen, der den Verlust gegentiber der Rate-Distortion- 
Schranke darstelit, welcher durch die suboptimalen kubischen 



l,53dB I Kaie-UislorUon-^cnran^e 
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1. StOrabstandsverlust bzgl. A » 6, 02 dB bei verschiedenen .4-Weiten 



Quantisierungszellen (gegentlber Kugeln ftlr D -4 oo) verur- 
sacht wird. In anderen Worten: Durch sphfirisch logarithmische 
Quantisierung ist es mttglich, den in Gl.(l) beschriebenen 
Verlust infolge der Kompandierung bis auf lediglich —1, 53 dB 
wieder zu kompensieren. Somit ist es durch die Wahl der Pa- 
rameters A und D theoretisch mOglich einen beliebig groBen 
Dynamikbereich (Bereich konstanten StOrabstandes Uber dem 
mittleren Signalpegel) zu erzielen, ohne nennenswerte Verlu- 
ste im maximal erzielbaren StOrabstand in Kauf nehmen zu 
mllssen. 

Abb. 2 zeigt als Beispiel den Abstand des StBrabstandes 
von der Rate-Distortion-Schranke R • 6dB abhangig vom 
mittleren Signalpegel (10 log 10 ( Varianz)) i.i.d. gauBverteilter 
Zufallsvariablen bei A = 48270 far unterschiedliche Dimen- 
sionality der spharisch logarithmischen Quantisierung. Diese 
imulationsergebnisse erfUUen genau die theoretische Analyse 
aB Abb. l.^arOberhlnaus wird deutlich, dass durch zwei 
kte mitVafchsender Dimensionalitat bei konstantem Para- 
Ster A zudem der Dynamikbereich stark ansteigt. Zum einen 
nimmt der Normierungsradius r 0 zu, zum anderen wirkt die 
Begrenzung der logarithmischen Kompression nur auf eine Di- 
mension, auf den Radius, wahrend fUr D — 1 Dimensionen die 
Proportionalitat der Ausdehnung der Quantisierungszelle zum 
Signalwert (=Bogenla*nge zu Radius bei fixer Winkeldifferehz) 
auch fUr kleinste Signalwerte erhalten bleibt. Somit erweitert 
sich der Dynamikbereich in D Dimensionen nSherungsweise 
auf 



B D « Bi + 20 log 10 (r 0 ) + 10 log 10 (^) . (28) 

Zudem steigt aufgrund des Mittelungseffektes innerhalb von 
jeweils D Werten die Obersteuerungsfestigkeit, wodurch eine 
weitere Vergr6Berung des Dynamikbereiches erzielt wird. 




mittlcrer Signalpegel [dB] 

Abb. 2. StCrabstand bei spharisch logarithmlBcher Quantisierung 
von i.i.d. gauOverteilten Signal werten 
{A = 48270, Bi - 93,7 dB, R > 1 bit/Abtaatwert) 



IV. KOMBINATION VON SPHARISCHER QUANTISIERUNG 
UND DPCM 

A. Problentstellung 

Die Differentielle Pulscodemodulation (DPCM) (siehe 
Abb. 3) Obertrfigt anstelle eines Abtastwertes q[k) selbst des- 
sen Differenz x[k] zu einem Prfidiktionswert Ziel ist die 
Erzeugung eines Differenzsignals nut kleinerer Varianz und 
geringerem Maximalwert als fUr die ursprCnglichen Abtast- 
werte i.d.R. um die gleiche Quantisierungsfehlerleistung wie 
bei PCM mit weniger Quantisierungsintervallen zu erreichen. 
Der PrSdiktionswert wird durch Filterung des zu quantisie- 
renden Signals mit einem Filter erzielt, durch das statistische 
AbhSngigkeiten bzw. eine 'Tarbigkeif des Quellensignals 
zur Redundanzreduktion ausgentitzt werden. Durch PrSdiktion 
lasst sich beispielsweise bei Audiosignalen oftmals bereits bei 
einer sehr kurzen FilterlSnge (und damit mit einem "wenig 
spezialisierten" Pradiktiorisfilter) ein Stflrabstandsgewinn von 
15dB und mehr erzielen. 

Es ist auBerdem anzumerken, dass durch DPCM keine 
systembedingte VerzOgerung des zeitdiskreten Signals ge- 
gentlber PCM hinzunehmen ist Weiterhin ist zu beachten, dass 
sich ein direkt messbarer Stbrabstandsgewinn (also ein Rein- 
gewinn bei gleicher Verzerrung des rekonstruierten Signals) 
bei DPCM nur in Form der sog. RUckwSrtspradiktion erzie- 
len IMsst; d.h. zur Gewinnung der Pradiktionswerte ist auch 
senderseitig die empfangsseitige Situation quantisierter Signal- 
werte $[k] anzuwenden, vgl. Abb. 3. Bei der anschlieBenden 
tibertragung ist zu beachten, dass der mittlere Signalpegel 
des Differenzsignals gegentlber dem Originalsignal um den 
Pradiktionsgewinn herabgesetzt ist (kleinere Varianz), sodass 
dieser Gewinn oftmals nicht vollstandig genutzt werden kann. 

Im Fall einer logarithmischen Quantisierung (SNR un- 
abhangig vom Signalpegel!) lassen sich diese beiden Gewinne 
jedoch tatsachlich addieren. An dieser Stelle tritt jedoch bei 




c 



Abb. 3. BlockKchaltbild: DPCM-Encoder mit ROckwtiitopriidlktion 



sphfirisch logarithmischer Quantisierung ein anderes Problem 
in den Vordergrund. Wie aus Abb. 3 deutlich wird, handelt 
es sich bei der RUckwartspradiktion urn ein rekursives Ver- 
fahren. Die quantisierten Werte q[k - i]\ i = 1,2, ... werden 
tch der Filterung vom darauffolgenden unquantisierten Wert 
subtrahiert und haben damit Einfluss auf die Quantisie- 
des unmittelbar nachfolgenden Wertes. Wahrend dies 
eindimensionalen Fall unproblematisch ist, steilt dies im 
mehrdimensionalen Fall und somit auch bei der spharischen 
Quantisierung einen grundsfitzlichen Widerspruch dar, weil 
hier jeweiis D Abtastwerte gemeinsam quantisiert und entspre- 
chend wieder gemeinsam rekonstruiert werden mtlssen. Bine 
gegenseitige Beeinflussung der Abtastwerte ist hier zun&chst 
nicht m&glich. 



B. Ldsungsgedanke 

Um den offensichtlichen Widerspruch der beiden Verfahren 
DPCM und mehrdimensionale Quantisierung aufzulttsen ist es 
erfordcrlich, fiir jeweiis D Abtastwerte diejenige Quantisie- 
rungszelle im D-dimensionalen spharisch quantisierten Raum 
zu finden, deren Rekonstruktionswerte q[k] im Sinne einer 
gemeinsamen Rekonstruktion von D Abtastwerten zu einen 
minimalen quadratischen Fehler (Quantisierungsfehler) in den 
Zeitschritten I mit k = I* D + * (i = 0,...,.D - 1) 




D-l 



d 2 («, q) = £ feP + 4[U ' I + i)) 2 



(29) 
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C Algorithnws 

Um die optimale Quantisierungszelle im D-dimensionalen 
spharischen Raum zu finden, verwenden wir einen iterativen 
Algorithmus, der - ausgehend von einem Startwert ("Start- 
zelle") - analog zu einem Gradientenabstiegsverfahren ar- 
beitet. Zur Ermittlung eines geeigneten Startwertes fiir die 
Abtastwerte q[D • J], ...,q[D • / + D - 1] wird von dem ak- 
tuellen Zustand des PrSdiktorfilters Hj>(z)z~ l zum Zeitschritt 
k = D • I eine Rekonstruktion der nSchsten D Werte gemaB 
einer Vorw&rtspradiktion ausgefiihrt, also bei Deaktivierung 
der Quantisierung (ADC und DAC in Abb. 3). Die auf diese 
Weise erhaltenen Werte x [D • I] , . . . , x[D • I + D - 1] am Eingang 
des Quantisierers werden nun .D-dimensional logarithmisch 
spharisch quantisiert und beschreiben damit die Startzelle 
(Zelle fUr 0-ten Iterationsschritt) fUr die Iteration. In jedem 



Iterationssschritt wird nun fUr alle Nachbarzellen der aktuellen 
Zelle bzw. fUr deren Rekonstruktionswerte $ [D • /), $[D • I + 
D-l] das Blockschaltbild durchlaufen, der Filterzustand von 
ifp(z)*- 1 wird dabei selbstverst&ndlich jedes Malzurttckge- 
setzt, und der jeweilige Quantisierungsfehler wird ermiltelt 
Hierbei ist zu bertlcksichtigen, dass beispielsweise fUr die 
Nachbarzellen in radialer Richtung i.d.R. andere Intervallan- 
zahlen fUr die Winkelquantisiemng vorliegen. Beispielsweise 
werden bei D » 3 die Azimutwinkel Qp\) bei groBen Eleva- 
tionswinkeln (<p2 nahe 90°) grtiber quantisiert als bei kleinen 
Blevationswinkeln nahe 0°), da bei letzteren Elevations- 
winkeln ein grOBerer Umfang des "Breitenkreises" auf der 
Kugeloberfl&che vorliegt. Diejenige Quantisierungszelle mit 
dem geringsten Quantisierungsfehler wird als "aktuell beste 
Zelle" ftlr den darauffolgenden Iterationsschritt verwendet, 

Simulationen zeigten, dass die im Mittel erforderliche Zahl 
an Iterationen bei ca. 0,25 (pro D Abtastwerte) liegt, und 
bereits bei einer BeschrSnkung auf max. 3 Iterationen keine 
signifikanten Verluste im Vergleich zu einem unbeschrHnkten . 
Suchraum zu verzeichnen sind, sodass die Suche nach der 
optimalen Quantisierungszelle fUr niedrige Dimensionszahlen 
durchaus in Echtzeit implementierbar ist Wie bereits aus 
Abb. 1 ersichtlich wurde, liegen die grCBten Gewinne ohnehin 
in diesem niedrigdimensionalen Bereich. 

Empfangsseitig ergibt sich keinerlei Anderung gegentiber 
Oblichen DPCM- Verfahren. An dieser Stelle soil darauf hin- 
gewiesen werden, dass die GesamtverzOgerung des Obertra- 
gungssystems lediglich D Abtastwerte betrSgt und somit ftir 
. Obertragungen, welche den Anforderungen einer extrem kur- 
zen VerzGgerung entsprechen mtlssen, hervorragend geeignet 
ist. 

D. Varianten des Verfahrens 

Zur Auffindung der gUnstigsten Quantisierungszelle die 
zu einer minimalen Verzerrung zJB. gemaB Gl. (29) fiihrt, 
sind alle Algorithmen der sog. Lattice-Decodierung bzw. zur 
Auffindung eines Maximum-Likelihood-Codeworts bei der 
Kanalcodierung prinzipiell zur Beschleunigung der Suche an- 
wendbar, vgl. z.B. [7] und die dort angegebenen Referenzen. 
Dabei ist vorzugsweise eine Transformation dieser Verfahren 
fiir Kugelkoordinaten vorzunehmen. 

Eine Variante ohne iterative Bestimmung der Quantisie- 
rungszelle ist dadurch gegeben, dass anstelle einer Berttck- 
sichtigung statistischer Abhangigkeiten innerhalb der aktu- 
ell zur Quantisierung anstehenden D Signalwerte durch ein 
lineares Pradiktionsfilter gemaB Abb. 3 (also DPCM), eine 
nicht gleichmafiige Quantisierung der Winkelkoordinaten vor- 
genommen wird. 

Dabei sind die linearen statistischen Bindungen der aktuell 
zur Quantisierung anstehenden D Signalwerte zu vergangenen 
Signalwerten mitteis klassischer DPCM mit Riickwartspradik- 
tion zu nutzen. Es entsteht also eine RUckwartspradiktion, bei 
der nach D Signalwerten eine Akttialisierung des Pradiktions- 
filters um D Schritte erfolgt. 

Es wird vorgeschlagen, die logarithmische Quantisierung 
des Betrags aufrecht zu erhalten, so dass anschlieBend die 



Signalpunkte so normiert wcrden kdnnen, dass sie auf ei- 
ner Kugel mit Radius 1 liegen. Dies erlaubt die anatyti- 
sche Berechnung der Wahrscheinlichkeitsdichtefunktion der 
Signalpunkte auf der Kugeloberfl&che, z.B. bei Annahme eines 
gauB'schen Signalprozesses aus der Autokorrelationsfunktion 
des Quellensignals, bzw. eine direkte experimentelle Ermitt- 
lung von relaliven H&ufigkeilen der Signalpunkte auf der 
Kugeloberflfiche. 

FQr M > 1 lasst sich hieraus z.B. mit Hilfe des Optimie- 
rungsansatzes nach [6] eine nichtgleichm&Bige kubische Quan- 
tisierung der Oberflfiche einer Kugel mit Radius 1 bestimmen. 
Dieser Ansatz besagt, dass im Mittel der Beitrag jeder Quan- 
tisierungszelle zum Quantisierungsger&usch gleich sein soil. 
Bezeichnet Ai die Kantenl&ige einer (D - l)-dimensionalen 
kubischen Quantisierungszelle sowie Zi den dazugehttrigen 
Mittelpunkt, also den empfangsseitigen Rekonstruktionsvektor 
dem Betrag 1, sowie f(tpu — i <Pd-i) die Wahrscheinlich- 
ichlefunktion bzw. relative HSufigkeit von Signalwerten 
der Kugeloberflfiche, so soil also nfiherungsweise erfUllt 
sein: 



A = 48270 bzw. A = 4858 und verschiedenen Dimensi- 
onszahlen D gemfiB Abb. 3 codiert und anschlieBend deco- 
diert. Der StBrabstand wurde jeweils durch Vergleich mit dem 
ursprilnglichen CD-Signal ennittelt (Abb. 4). Da die hier ein 
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Hieraus folgt unmittelbar eine (D — l)-dimensionale 
Kompressor-Funktion fc(y>i,...,y>D-i). durch die die 
ungleichmSBige Quantisierung der Kugeloberflfiche eipdeutig 
bestimmt wird. Die ungleichmfiBige Quantisierung kann 
dann beispielsweise wie im eindimensionalen Fall durch eine 
nichtlineare Verzerrung des Vektors *=(y>i, <PD-i) in einen 
Vektor z'=k(z) t anschlieBende gleichmfiBige Quantisierung 
zum Vektor z\ gemfiB Abschnitt HI-C und Anwendung der 
Umkehrfunktion nachfolgend zur Gewinnung 

des Rekonstruktionsvektors vollzogen werden. 

r orzugsweise ist die resultierende (D — l)-dimensionale 
pressorfunktion durch eine analytisch beschreibbare 
:ion anzunahern, oder, analog zu der bei der eindi- 
mensionalen Kompression Ublichen Darstellung mittels Ge- 
radenstQcke (z.B. 13-Segment-Kennlinie, siehe[3]) durch 
(D — 2)-dimensionale Teilebenen (mit konstanten partiellen 
Ableitungcn) zu approximieren. Bei einer iiber dem Raum 
(-f,!]^" 2 x (-7r,7r] orthogonalen Rasterung dieser Ap- 
proximation Uber den D — 1 Winkeln kann eine einfache 
Auswertbarkeit der Kompressorfunktion erreicht werden. 

V. SlMULATIONSERGEBNISSB 

Zur Simulation von sphfirisch logarithmischer Quantisierung 
in Kombination mit DPCM gemfiB Abschnitt IV-C wurde die 
Quantisierung der Arie "Der Vogelffinger bin ich ja" von Wolf- 
gang Amadeus Mozart [5] untersucht Es weist einen mittleren 
Signalpegel von -32,15dB auf. Diese Arie bietet aufgrund 
von Signaldynamik und Klangfarben (Vorspiel, Gesang, Rohr- 
flete) ein anspruchsvolles Beispiel zur Audiosignalcodierung. 
Es wurde fUr R = 4 und R = 7 Bit pro Abtastwert bei 




R = 7 bit/ Abtastwert 
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Abb. 4. SNR bei sphflrisch logarithnuscher Quantisierung 
in Kombination mit DPCM 
M Der Vogelftnger bin ich ja" (W. A. Mozart). 
Spielzeit « 3:02 min, -32, 15 dB , P - 2 
R = 7 bit/ Abtastwert , A = 4858 
R = 4 bit/ Abtastwert , A = 48270 

universeller Prfidiktor von geringer Prftdiktionsordnung (P = 
2) verwendet wurde, kttnnen diese Simulationsergebnisse als 
stellvertretend ftlr eine Vielzahl von Audiosignalen betrachtet 
werden. Abb. 4 zeigt die gemessenen StOrabsta*nde, gemittelt 
Uber das ganze MusikstUck. Bei geringen Raten ergeben sich 
aufgrund einer verbesserten Mittelwertbildung bei steigender 
Dimensionszahl D sogar weit grBBere StBrabstandsgewinne 
als nach Abb. 1. Beim Vergleich der Ergebnisse sind die 
unterschiedlichen Werte fUr den Parameter A zu beachten. 

VI. ZUSAMMENFASSUNG 

Es wurde ein wellenformerhaltendesDigitalisierungsverfah- 
ren fiir analoge Quellensignale vorgestellt, das zum einen 
den Gewinn durch mehrdimensionale Quantisierung mit den 
Vorteilen der logarithmischen Quantisierung verkniipft und 
zum anderen in der Lage ist den objektiven Stttrabstand durch 
Addition von Prfidiktionsgewinnen zusStzlich zu erhOhen. Ab- 
schlieBend sei angemerkt, dass sich dieses Verfahren neben 
einem gtinstigen Austausch von Rate und Verzerrung vor allem 
durch einen extrem grofien Dynamikbereich und dabei einer 
BuBerst geringen systembedingten Verztigerung des Signals um 
wenige Abtastintervalle auszeichnet. 
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AnsorQche 



Verfahren zur Digitalisierung analoger Quellensignale, mit dem Schritt: 
- DurchfQhren einer D dimensionalen logarithmischen spharischen Quan- 
tisierung. 




2. Verfahren zur differentiellen Pulscodemodulation, 

DurchfQhren einer D dimensionalen logarithmischen spharischen Quan- 
tisierung. 



3. Verfahren nach Anspruch 1 oder 2, wobei 

eine logarithmische Quantisierung des Betrages erfolgt. 



Verfahren nach Anspruch 2, mit den Schritten 

DurchfQhren einer Forwartspradiktion zur Ermittlung eines Startwertes 
fQr Abtastwerte der Quantisierung auf der Basis des aktuellen Zustandes 
eines Pradikatorfilters; • 

Rekonstruktion der nachfolgenden D Abtastwerte, 



-2- 



DurchfQhren einer D dimenslonalen logarlthmlschen spharischen Quan- 
tisierung der durch die Forwartspradiktion erhaltenen Werte, urn eine 
Startzelle zu bestimmen, 

iteratives Durchlaufen der Pradiktion der differenziellen Pulscodemodu- 
lation, urn eine Quahtisierungszelle mit dem geringsten Quantisierungs- 
fehler zu bestimmen. 

5. Verfahren zur Digitalislerung analoger Quellensignale, insbesondere nach 
einem der vorherigen AnsprQche, wobei 

eine Lattice-Decodierung zur Bestimmung einer gOnstigen Quantislerungs- 
zelle durchgefQhrt wird. 

6. Verfahren nach Anspruch 5, wobei 

eine Transformation in Kugelkoordinaten durchgefQhrt wird. 

7. Verfahren zur Digitalisierung analoger Quellensignale, insbesondere nach 
einem der vorherigen AnsprQche, wobei 

eine ungleichmafiige Quantisierung der Winkelkoordinaten der sphari- 
schen oder Kugelkoordinaten durchgefQhrt wird. 

8. Verfahren nach Anspruch 7, wobei 

eine RQckwartspradiktion einer differentiellen Pulscodemodulation 
durchgefQhrt wird, wobei nach der Verarbeitung von D Abtastwerten ei- 
ne Aktualisierung des Pradiktionsfilters urn D Schritte erfolgt. 

9. Verfahren nach Anspruch 7 oder 8, wobei 

der Betrag logarithmisch quantisiert wird. 

1 0. Vorrichtung zur Digitalisierung analoger Quellensignale, mit 

- einer D dimensionalen logarithmischen spharischen Quantisierungsein- 
richtung. 

.11. Vorrichtung zur differentiellen Pulscodemodulation, 
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elner D dimenslonalen logarithmischen spharischen Quantislerungseln- 
rlchtung. 

1 2. Vorrichtung nach Anspruch 1 0 oder 1 1 , wobei 

elne logarlthmlsche Quantislerung des Betrages erfolgt. 

1 3. Vorrichtung nach Anspruch 1 1 mit, 

einer Forwartspradiktionselnrichtung zur Ermittlung eines Startwertes fur 
die Abtastwerte der Quantislerung auf der Basis des aktuellen Zustan- 
des eines Pradikatorfilters, 

einer Rekonstruktionseinrichtung zur Rekonstruktion der nachfolgenden 
D Abtastwerte, 

einer D dimenslonalen logarithmlschen spharischen Quantisierungsein- 
richtung zur Quantislerung der durch die Forwartspradiktion erhaltenen 
Werte, um eine Startzelle zu bestimmen, 

wobei die Pradiktion der differenzlellen Pulscodemodulation iterativ 
durchiaufen wird, um eine Quantisierungszelle mit dem geringsten 
Quantisierungsfehler zu bestimmen. 

1 4. Vorrichtung zur Digitalisierung analoger Quellensignale, insbesondere nach 
einem der vorherigen AnsprQche, wobei 

eine Lattice-Decodlerungseinrichtung zur Bestimmung einer gdnstigen 
Quantisierungszelle durchgefuhrt wird. 

15. Vorrichtung nach Anspruch 14, mit 

einer Transformationseinrichtung zur Durchfuhrung einer Transformation 

in Kugelkoordinaten. 

1 6. Vorrichtung zur Digitalisierung analoger Quellensignale, insbesondere nach 
einem der vorherigen AnsprQche, wobei 

eine ungleichmaliige Quantisierung der Winkelkoordinaten der sphari- 
schen Oder Kugelkoordinaten durchgefuhrt wird. 

17. Vorrichtung nach Anspruch 16, wobei 



eine RQckwartspradiktion einer differentiellen Pulscodemodulation 
durchgefOhrt wird t wobei nach der Verarbeitung von D Abtastwerten ei- 
ne Aktuallsferung des Pradiktionsfilters urn D Schritte erfolgt. 

Vorrichtung nach Anspruch 16 oder 17, wobei 
der Betrag logarithmisch quantisiert wird. 



